In silico screening of some compounds derived from the desert medicinal plant Rhazya stricta for the potential treatment of COVID-19.

The latest coronavirus pandemic (SARS-CoV-2) poses an exceptional threat to human health and society worldwide. The coronavirus (SARS-CoV-2) spike (S) protein, which is required for viral-host cell penetration, might be considered a promising and suitable target for treatment. In this study, we utilized the nonalkaloid fraction of the medicinal plant Rhazya stricta to computationally investigate its antiviral activity against SARS-CoV-2. Molecular docking and molecular dynamics simulations were the main tools used to examine the binding interactions of the compounds isolated by HPLC analysis. Ceftazidime was utilized as a reference control, which showed high potency against the SARS-CoV-2 receptor binding domain (RBD) in an in vitro study. The five compounds (CID:1, CID:2, CID:3, CID:4, and CID:5) exhibited remarkable binding affinities (CID:1, - 8.9; CID:2, - 8.7; and CID:3, 4, and 5, - 8.5 kcal/mol) compared to the control compound (- 6.2 kcal/mol). MD simulations over a period of 200 ns further corroborated that certain interactions occurred with the five compounds and the nonalkaloidal compounds retained their positions within the RBD active site. CID:2, CID:4, and CID:5 demonstrated high stability and less variance, while CID:1 and CID:3 were less stable than ceftazidime. The average number of hydrogen bonds formed per timeframe by CID:1, CID:2, CID:3, and CID:5 (0.914, 0.451, 1.566, and 1.755, respectively) were greater than that formed by ceftazidime (0.317). The total binding free energy calculations revealed that the five compounds interacted more strongly within RBD residues (CID:1 = - 68.8, CID:2 = - 71.6, CID:3 = - 74.9, CID:4 = - 75.4, CID:5 = - 60.9 kJ/mol) than ceftazidime (- 34.5 kJ/mol). The drug-like properties of the selected compounds were relatively similar to those of ceftazidime, and the toxicity predictions categorized these compounds into less toxic classes. Structural similarity and functional group analyses suggested that the presence of more H-acceptor atoms, electronegative atoms, acidic oxygen groups, and nitrogen atoms in amide or aromatic groups were common among the compounds with the lowest binding affinities. In conclusion, this in silico work predicts for the first time the potential of using five R. stricta nonalkaloid compounds as a treatment strategy to control SARS-CoV-2 viral entry.

Pharmacophore-Based Virtual Screening, Quantum Mechanics Calculations, and Molecular Dynamics Simulation Approaches Identified Potential Natural Antiviral Drug Candidates against MERS-CoV S1-NTD.

Middle East respiratory syndrome coronavirus (MERS-CoV) is a highly infectious zoonotic virus first reported into the human population in September 2012 on the Arabian Peninsula. The virus causes severe and often lethal respiratory illness in humans with an unusually high fatality rate. The N-terminal domain (NTD) of receptor-binding S1 subunit of coronavirus spike (S) proteins can recognize a variety of host protein and mediates entry into human host cells. Blocking the entry by targeting the S1-NTD of the virus can facilitate the development of effective antiviral drug candidates against the pathogen. Therefore, the study has been designed to identify effective antiviral drug candidates against the MERS-CoV by targeting S1-NTD. Initially, a structure-based pharmacophore model (SBPM) to the active site (AS) cavity of the S1-NTD has been generated, followed by pharmacophore-based virtual screening of 11,295 natural compounds. Hits generated through the pharmacophore-based virtual screening have re-ranked by molecular docking and further evaluated through the ADMET properties. The compounds with the best ADME and toxicity properties have been retrieved, and a quantum mechanical (QM) based density-functional theory (DFT) has been performed to optimize the geometry of the selected compounds. Three optimized natural compounds, namely Taiwanhomoflavone B (Amb23604132), 2,3-Dihydrohinokiflavone (Amb23604659), and Sophoricoside (Amb1153724), have exhibited substantial docking energy >-9.00 kcal/mol, where analysis of frontier molecular orbital (FMO) theory found the low chemical reactivity correspondence to the bioactivity of the compounds. Molecular dynamics (MD) simulation confirmed the stability of the selected natural compound to the binding site of the protein. Additionally, molecular mechanics generalized born surface area (MM/GBSA) predicted the good value of binding free energies (ΔG bind) of the compounds to the desired protein. Convincingly, all the results support the potentiality of the selected compounds as natural antiviral candidates against the MERS-CoV S1-NTD.